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Abstract. We examine the effects of cosmic strings on structure formation and on the 
ionization history of the universe. While Gaussian perturbations from inflation are known 
to provide the dominant contribution to the large scale structure of the universe, density 
perturbations due to strings are highly non-Gaussian and can produce nonlinear structures 
at very early times. This could lead to early star formation and reionization of the universe. 
We improve on earlier studies of these effects by accounting for high loop velocities and for the 
filamentary shape of the resulting halos. We find that for string energy scales Gfi > 10^^, the 
effect of strings on the CMB temperature and polarization power spectra can be significant 
and is likely to be detectable by the Planck satellite. We mention shortcomings of the 
standard cosmological model of galaxy formation which may be remedied with the addition 
of cosmic strings, and comment on other possible observational implications of early structure 
formation by strings. 
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1 Introduction 

Cosmic strings are linear topological defects that could be formed at a phase transition in 
the early universe [1] . They are predicted in a wide class of particle physics models. Some 
superstring-inspired models suggest that fundamental strings may also have astronomical 
dimensions and play the role of cosmic strings [2-4] . 

Strings can be detected through a variety of observational effects. Oscillating loops of 
string emit gravitational waves - both bursts and a stochastic background. They can also 
be sources of ultrahigh-energy cosmic rays. Long strings can act as gravitational lenses [5- 
7] and can produce characteristic signatures in the CMB: discontinuous temperature jumps 
[8, 9] along the strings (Gott-Kaiser-Stebbins effect) and a B-mode polarization pattern 
due to rotational perturbations induced by strings [10-12]. Formation and evolution of 
cosmic strings and their possible observational effects have been extensively discussed in the 
literature; for a review and further references see, e.g., [13-15]. 

The strength of gravitational interactions of strings is characterized by the dimensionless 
number Gfx, where // is the mass per unit length of string and G is Newton's constant. Early 
work on cosmic strings was largely motivated by the idea that oscillating string loops and 
wakes formed by rapidly moving long strings could serve as seeds for structure formation 
[16-18]. This scenario, which required Gfi ~ 10~^, has been conclusively ruled out by CMB 
observations. Present CMB bounds constrain G/j, to be less than about 2 x 10~'' [19-24]. A 
much stronger bound, Gfi < 4 x 10~^, has recently been claimed in [25] (but see also [26]) 
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and [27] based on the millisecond pulsar observations of the gravitational wave background. 
(We shall argue in section 7 that this bound could be significantly relaxed, in the light of the 
latest string simulations.) 

The CMB and other data are consistent with a Gaussian spectrum of density perturba- 
tions predicted by the theory of inflation; the contribution of strings, if any, can account for 
no more than 2% of the spectral power. This does not mean, however, that strings always 
played a subdominant role in structure formation. Density perturbations due to strings are 
highly non-Gaussian and can produce nonlinear structures at very early times. This could 
result in early star formation and reionization of the universe [28-33]. Such early structure 
formation could manifest itself in CMB temperature and polarization [34] spectra and could 
also produce a detectable 21-cm signal [35-38]. 

These and other effects of strings depend on the details of string evolution, which until 
recently remained rather uncertain. Early string simulations suggested that loops produced 
by the string network are short-lived and/or very small^ [39-43]. Hence, it was assumed 
in Refs. [28-31, 34-37] that the main effect on structure formation was due to long string 
wakes. However more recent high-resolution simulations have demonstrated that a substan- 
tial^ fraction of the network energy goes into large non-self-intersecting loops, having length 
of about 5% of the horizon [44-46]. This pattern is established only after a long transient 
period dominated by very small loops; apparently it is this period that was observed in pre- 
vious simulations. Ref. [32] studied structure formation and reionization in this new string 
evolution scenario, but some important aspects of that scenario had not yet been recognized 
at the time. 

In the present paper we shall reexamine early structure formation by cosmic strings and 
its effect on the ionization history of the universe. An important fact not taken into account 
in Ref. [32] is that string loops are typically produced with high velocities: v ~ 0.3 for the 
largest loops and even higher for smaller ones [46]. Accretion on such rapidly moving loops 
is rather different from spherical accretion on a stationary point mass which was assumed in 
[32] . We also make use of the latest string simulations which yielded more reliable results for 
the size and velocity distributions of cosmic string loops. 

The paper is organized as follows. In the next section we review the evolution of cosmic 
strings and introduce the relevant loop distributions. In section 3 we study accretion of 
matter onto a moving loop, first treating the loop as a point mass and then accounting for 
its finite size (which turns out to be significant). We find that loop seeded halos have the 
form of highly elongated filaments which then fragment into smaller "beads" . The resulting 
halo spectrum is calculated in section 4. The baryon collapse fraction in string-seeded halos 
is calculated in section 5 and is then used in section 6 to study the reionization history of 
the universe and the possible effect of strings on the CMB temperature and polarization 
spectra. Finally, in section 7 we comment on possible observational implications of the early 
structure formation by cosmic strings. In appendix A we discuss the "rocket effect" - the 
self-acceleration of string loops due to asymmetric emission of gravitational waves - and show 
that it is not significant for structure formation. 



^The simulation reported in [39] found a scaling loop population including large but fragmenting loops. 
Subsequent analytic work [40] assumed small loop production to be more significant for the loop energy 
density. 

^The discrepancy between this result and the results of [40] will be addressed in a future publication. 
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2 Cosmic string evolution 



2.1 Long strings 

An evolving network of cosmic strings consists of two components: long strings and sub- 
horizon closed loops. The long string component evolves in a scaling regime in which the 
typical distance between the strings and all other linear measures of the network remain at 
a constant fraction of the horizon size (ih(i)- This dynamic is maintained by the produc- 
tion of loops via reconnection at string intersections, and subsequent evaporation of loops 
by gravitational radiation. On super-horizon length scales, long strings have the form of 
random walks, while on smaller scales they exhibit a wide spectrum of wiggles and kinks, 
which are remnants of earlier reconnections down to the scale of initial formation. These 
small-scale features are gradually smoothed out by expansion, small loop production, and by 
gravitational radiation back-reaction. 

The total energy density in the long ("infinite") strings can be expressed as 



where /i is the mass per unit length of string (equal to the string tension), the radiation 
era horizon distance is = 2t, 7 is a constant coefficient, and the quantity jd^ defines 
the average inter-string distance. The idea that the network properties are determined by 
the horizon distance alone is known as the scaling hypothesis [1, 47]; for the long string 
component it has been confirmed through numerical simulations [41, 42, 48]. 

The coefficient 7 in (2.1) depends on the reconnection probability prec of intersecting 
strings. For "ordinary" particle physics strings, Prcc = and simulations give 7 ~ 0.15 in 
the radiation era. For cosmic superstrings, the reconnection probability is expected to be 
much smaller, 0.1 ^ Prec 

> 10-3 [49], resulting in a smaller value of 7 and a denser string 

1/2 

network. Simple arguments suggest that 7 cx pric [50, 51], but simulations indicate a weaker 
dependence [52]. Here we shall focus on ordinary strings'^ with prec = 1- 

2.2 Loop distribution: a simple model 

The length distribution of loops produced by the network (the so-called loop production 
function) obtained in recent simulations has a double-peak structure, indicating two different 
populations of loops. First, there is a scaling loop distribution, with a typical loop length 
of about 5% of the horizon and a wide tail extending to smaller scales. Then there is a 
non-scaling distribution of very small loops with sizes comparable to the initial scale of the 
network at formation. The hight of the non-scaling peak is observed to decrease somewhat 
in the course of the simulation, with the extra power contributing to the short-length tail 
of the scaling peak. One might expect that the non-scaling peak will eventually disappear 
[44-46], but Ref. [56] argued that it could also survive at late times. In the latter case, 
the asymptotic double-peak loop distribution will still scale, but the typical loop size in the 
short-length peak will be set by the gravitational radiation damping. 

At present one can only guess which of the two options will be supported by future 
simulations, but fortunately this is not important for the purposes of the present paper. The 

^We assume that strings are accurately described by the Nambu-Goto equations of motion. Abehan field 
theory simulations (see [53] and references therein) yield different network properties, suggesting that the main 
energy loss mechanism of the network is direct production of massive particles from long strings. However, 
analysis in Refs. [54, 55] indicates that massive particle production is a transient phenomenon, caused by 
short-wavelength excitations in the initial state. 
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non-scaling loops are highly relativistic, with the dominant part of their energy being kinetic 
energy. This energy redshifts with the expansion and has very little effect on the gravitational 
clustering. The same applies to loops in the tail of the scaling distribution. Only loops near 
the peak of the scaling distribution, which are formed with mildly relativistic velocities, are 
relevant for structure formation. Moreover, halos formed by small loops at the tail of the 
distribution have small virial velocities. This inhibits cooling and star formation in such 
halos. 

The string loops of interest to us here will be those which formed during the radiation 
era but have not yet decayed at teq- The energy density of loops that were chopped off the 
network in one Hubble time is comparable to the energy density of long strings. However, 
the loop energy redshifts like matter, while the long string energy redshifts like radiation. 
So, if loops are long and live much longer than a Hubble time, they dominate the energy of 
the network and play dominant role in structure formation. 

Using energy conservation, the power flowing into loops per unit physical volume during 
the radiation era must obey 

F— sloops — , Hoc — ,3 ) {■^■■^J 

t 

where (f^) is the mean square velocity of the long strings. Simulations give (u^) ~ 0.4 and 
V := /i~^d^/9_5.ioops ~ 50 in this epoch. Neglecting (for the time being) the center-of-mass 
loop velocity, the comoving number density distribution of loops n{t,m) should satisfy 



1 d 



a^{t) dt \Jq 



00 



n{t,m)mdm\ = p^loops- (2.3) 



Here a{t) is the scale factor, and the comoving loop distribution is related to the physical 
distribution by n{t,m) = a'^(t)nphys(i, ^T^)• 

Since we are only interested in loops near the peak of scaling distribution, it is a rea- 
sonable approximation to assume that all relevant loops are formed having length equal to a 
fixed fraction of the horizon size. Hence, we set the loop mass at formation to be 

m = afj,d\i, (2-4) 

with a ~ 0.05. The distribution of such loops obeys 

'n{t,m)m] = 5{m — 2fiat)-^-!^, (2-5) 



a3(t) dt ^ J ^ ^ ^ 8t3 ' 

and so 

and n{t, m) = otherwise. Here, 6V ^ 7 reflects the fact that the effective power flowing into 
these large loops is substantially less than the total power P 51 suggested by simulations 
[46]. We have also used a{t) = {t/t^qY^'^, normalizing to a{teq) = 1- 

The apparent divergence of loop mass density mn{t,m)dm at small m is absent if we 
include the decay of loops due to gravitational radiation, m ~ —TGfj?, with T ~ 50. This 
correction is only significant for very small loops. 
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2.3 Incorporating loop speed 

We now perform a more precise calculatfon which includes finite velocity effects, which are 
significant. We characterize the number density of loops by their mass m and speed v. All 
quantities will be expressed per comoving volume from here onward. The comoving number 
density of cosmic string loops with mass between m and m + dm and with speed between v 
and v + dv is given by n{t, m, v)dv dm. The loop production function is g{t, m, v)dt dv dm., the 
number of such loops per comoving volume produced by the cosmic string network between 
time t and t+dt. The redshifting of loop speeds is given by ~ —Hv where H = a/ a = l/(2t) 
is the Hubble rate, and we are assuming non-relativistic f , which is justified for the subset 
of loops which contribute to star formation. 

We can integrate the production rate to find the number density per comoving volume 

dvn! dv' 

nit,m,v) = dt'g{t\m\v') — — . (2.7) 

This states that the number of loops of mass m and speed v at time t is the integral of the 
production rate of loops over all prior times t' < t, where the relevant production is of loops 
which will eventually have mass m and speed v. Hence the first step is to write down the 
solution to the flow, namely 

m'{t'; m,t) ^ m + TGi?{t - t'), (2.8) 
v'it';v,t)^v^ (2.9) 

These tell us e.g., what speed v' a loop must have at production time t' in order for it to 
have a speed v at time t. The Jacobian factor ^ captures the changing size of the volume 
element dv. The same applies to m as well. We are neglecting the effects on v of anisotropy 
in the loop gravitational radiation, also known as the "rocket effect" , which causes the loops 
to accelerate as they evaporate. This will be justified in the Appendix. 
The energy of a loop is given by m/Vl — v'^, and so using (2.2), 

r/" d„.,„. (2.10) 

<(«) Jo Jo aHt) 

Since all quantities of interest are linear in g, we can maintain generality while treating g as 
a delta-function source, i.e., 

a^(t) 8Pii^/\ — f 2 

g{t, m, v) = ^3^^^^ 5 (m - aiid\,{t)) 5 {v - a^) , (2.11) 

where in this form we can imagine taking the values suggested by simulations, namely 

a = 0.05, = 0.3, 5V = 7. (2.12) 

The results using the actual loop production function, including small relativistic loops, 
can always be obtained by integrating any expression involving SV with the replacement 
— 7'(numerical loop-production distribution) using simulation data from e.g., [46]. Because 
of this substitution possibility, we will refer to (2.11) as the delta-function form of g, and 
(2.12) as the delta- function approximation for g. 
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Combining the equations of this section, the comoving loop number density is 



n{t,m,v) = I dt'g{t',m,v')-^ — , (2.13) 
* , .a?{t')6Vfi , , , , , „^ , / , ^ a(t) 



3 Accretion onto a cosmic string loop 

As we explained in Section 2.2, the dominant loop energy density will always be from loops 
produced in the radiation era. This is because the integrated power flowing into loops scales 
like 1/t^, whereas the subsequent dilution only scales like 1/a^. Thus in the radiation-era 
(when a ~ t^^"^), loops pile up from early times. This growth is cut off by loop evaporation 
to gravitational radiation, which is a slow process for low tension strings. So structure 
formation is mainly sensitive to radiation-era loops which have survived to the time of matter 
domination. 

Of particular interest are halos which become large enough for stars to form [57]. This 
requires the baryons to collapse to sufficient density, which can only happen due to dissipation. 
When the neutral hydrogen atoms have large enough virial velocity, collisions will have 
sufficient center-of-mass energy to excite their electrons and radiate. The energy escaping 
with the radiated photons ensures that the hydrogen loses gravitational potential energy. This 
cooling moves the hydrogen toward the center of the halo. The critical virial temperature 
for efficient cooling can be as low as 200K in the case of molecular hydrogen, but we will 
neglect this mechanism since the molecule is too fragile to survive UV light from the very 
first stars. Instead, we will consider two possible values, lO^K and lO^K, below which star 
formation does not occur. ^ Because the virial temperature is proportional to a positive power 
of the halo mass, which is proportional to the seed (loop) mass, there is a minimum loop 
size which can lead to star formation. We can conclude that early stars can only come from 
loops produced late enough to have this size. We will denote the minimum scale factor after 
which star-forming loops can be produced by a*, which we will now estimate. Throughout, 
we define the scale factor at matter-radiation equality to be unity, a(teq) = 1. 

3.1 Spherical accretion 

For a spherical halo of mass M formed at redshift z, the virial temperature is given by [57] 
and so the minimal mass of a halo that has T^\^ > 10^ K at redshift z is 

3/2 



M.(z) = 3xlO^Mo(l + z)-^/^(^y . (3.2) 



^Because of the relative motion between the baryon and dark matter fluid, smaller halos will not lead to 
significant star formation [58, 59]. 
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Because the vast majority of loops are moving rather quickly, the more appropriate 
bound on production time considers elongated filaments, rather than spherical halos. Elon- 
gated halo filaments have a virial temperature [60] dependent entirely on their linear mass 
density /ifii, which we will calculate below. These filaments will later collapse into beads, 
which subsequently merge into larger beads, but this process does not significantly affect the 
virial temperature. 

3.2 Accretion onto a moving point mass 

If we assume the loop is non-relativistic and has velocity feq > in the +?/-direction at teq, 
its subsequent velocity is given by 

via) = ^, (3.3) 
a 

where a = (t/teq)^^'^- The trajectory in comoving coordinates is then 

y{a) = Sveqt^q (^1 - -^^ . (3.4) 
Using our delta-function loop production function, the loop velocity at equality is 

t^eq ~ a^tti, (3.5) 

where Oj < 1 is the scale factor at loop production, relative to a^q = 1. 

Given a point mass in the above trajectory, we can find the cylindrically symmetric 
turn-around surface by considering a particle at comoving initial location {xo,yo,0). (Here 
we closely follow the analysis in [61].) Let us label the moment of closest approach of the 
mass by the scale factor ao, i.e., the loop trajectory y(a) obeys y(ao) = yo- Using the impulse 
approximation, the velocity kick on the particle due to the passing point mass is 

Gm Gm 

Vx ~ 7 Tn Ato ' (3.6) 

[aoxoy aoXov{ao) 

where 

Ato~^. (3.7) 
v{ao) 

— 1 3/2 

This particle will then be displaced after one subsequent Hubble time Hq ~ icq^o 
amount 

aoAj;o ~ VxieqOo^^' (3-8) 

and so the corresponding density perturbation is 

Axo Gmtcq , . 

6o^- 1/2 ' (3-9) 



^0 x^v{ao)aQ 



which grows to be 



The turnaround surface is given by 6 {a) ~ 1, so the profile of the collapsed region from a 
passing point mass with arbitrary velocity v{ao) is 

a;ta(ao,a) 3/2 ^ (3.11) 

where the time dependence is in a = (t/tgq)^^^, and the y-dependence is in ao(yo) via (3.4) 
with y{ao) = hq. Notice a > gq, since the turnaround surface extends only behind the loop's 
^/-position. The halo mass is then given by the total mass inside the turnaround surface, 

M{t) = pa\t) r^'''^'^\xl{aoiyo),a{t))dyo = pa\t) f ^xlfi^^^^dt, (3.12) 







a.,N 7i^4Xao) 1 /"* TrGmteqoW rna{t) /"* tl^^ 

pa [t) / dto = ^_„,2 / = -^r- / ^^*o 



an 67rGt2 / 5/2 a I ,5/3 

= J Ht) - 1] . 

Interestingly, the loop trajectory does not affect the halo mass at this level of approximation. 
A more accurate analysis by Bertschinger [61], using the Zel'dovich approximation, gives 
M{t) = ^ma{t). We will simply use 

M(t)~ma(t). (3.13) 

We can solve for the shape and size of the turnaround surface by combining (3.3), (3.4) 
& (3.11). This gives 



Xta.{y,a) = teqW ( ^ ~ o ^~r~] < y < Su^aiteq ( I ^ ) . (3.14) 

Notice that these halos are very elongated: the eccentricity after a Hubble time is given 

by 



a~2 



3/2 1/2 

~ lOVlf . (3.15) 

^a(jp 



Here we have assumed ai = ^Jtij tcq ~ 0.1, and used the delta- function approximation 
~ 0.3 and a ~ 0.05. We refer to these elongated structures as filaments. Halos will form 
from linear instabilities of the filaments. 

In this discussion we disregarded the so-called rocket effect - the self-acceleration of the 
loop due to asymmetric emission of gravitational waves. This has negligible effect on the 
loop's velocity, except toward the end of the loop's life, when the loop can be accelerated to 
a mildly relativistic speed, v < 0.1. The loop trajectory with the rocket effect included is 
discussed in the Appendix, where it is shown that the rocket acceleration has little influence 
on halo formation. 

3.3 Accretion onto a finitely extended loop 

Since the loop has a finite radius R = f3m/p, with f3 ~ 0.1, only the matter outside this 
distance will feel a momentum kick from the passing loop.^ We should then only consider 



^The rapidly oscillating string will leave wakes of overdensity behind the fast segments, i.e., even inside 
the loop radius R. We neglect this effect, since only a small fraction ~ SttGh/v of the material is affected, 
and the wakes are probably too thin to allow star formation. 
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the portion of the turnaround surface where 



(3.16) 



or using (3.3) & (3.11), 



< 



ao 



2 ' 

i 



(3.17) 



where again, oq should be thought of as a measure of the y-coordinate given by yo = y{ao) 
above. This equation is a restriction on the vahdity of the turnaround surface (3.14); in 
places where the turnaround surface is smaller than the comoving loop radius, it does not 
exist and so should be thought of as ending rather than closing. 

For the faster loops (whose turnaround surfaces have smaller physical radius c.f. (3.11)), 
the loop radius will entirely cloak the turnaround surface at the onset of matter domination. 
Because the turnaround surface radius grows relative to the comoving loop radius, there will 
eventually be a time when the surface emerges. The surface will emerge as a hoop surrounding 
the loop, which grows into a cylinder. The front of the cylinder will move forward with the 
loop, and continue to grow in diameter with expansion. The back of the cylinder will extend 
back to the location where the turnaround surface is smaller than the comoving loop radius. 
Eventually the back of the cylinder will reach y = 0, i.e., the location of the loop at the onset 
of matter domination. 

We can now talk of three distinct phases of loop-seeded filaments in the matter era. The 
early type are those which are not accreting, because no part of the turnaround surface has 
emerged from behind the loop radius. The late type are filaments whose entire turnaround 
surface is larger than the loop radius. Because the point-mass approximation holds in this 
case, the total filament mass is given by (3.13). We will call these "normal growth" filaments. 
The intermediate type of filament are those whose turnaround surfaces are growing both in 
the positive y-direction with the loop motion, as well as in the negative y-direction as more 
and more of the turnaround surface emerges from beneath the loop radius. We will call 
these "accelerated growth" filaments, since they are catching up from having zero mass to 
eventually have normal mass. 

We can find when the turnaround surface first emerges from behind the loop radius R 
by setting oq = a in (3.17) to find 



a > On 



Gil 



2/5 



(3.18) 



The turnaround surface stops its accelerated growth when (3.17) is satisfied all the way back 
to yo = 0, i.e. oq = 1, and so the normal-growth regime for the filament takes over after 



5/2 



(3.19) 



The turnaround surfaces for a = 10, 20. ..370 are shown in figure 1. 

Given these turnaround surfaces, the total filament mass can be shown to be 



ry{t) 

Mm{t) = pcq / x{t)^Jyo 



ro 



. 5/3 



mm 

m(a — 1) 



Cl < Omin 
O-min < ffl < O-n 

O-max ^ 



(3.20) 
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Comoving turnaround surfaces and loop radius 



l.Ot 



0.8 i 



0.6 1 




0.0000 



0.0005 



0.0010 0.0015 
/ ^eq 



0.0020 



0.0025 



Figure 1. Time-lapse (Aa = 20) illustration of comoving turnaround surfaces (thin blue lines) for a 
loop of tension G/i = 10^^. The comoving loop extension is shown dotted in red. Notice the aspect 



ratio is reduced by a factor of 400:1, and units are <oq- One can read off a^^^ ~ 10, Oi, 



370 from 



the figure. We use at = 1/2 and the delta-function approximation parameters of (2.12) to determine 
the mass and speed of the loop. The rocket effect is included for completeness (but see Appendix). 



Note that in the delta-function approximation for g where oc oc m, the filament mass 
in the accelerated growth regime is the same for all loops, independently of their mass m. 
With the loop trajectory from (3.4), the physical length of the filament is given by 



r 

3a a t (^ 
^ Sa^aiteqia - ^/a) 



C^min ^ CL <C Q>max 



(3.21) 



a 



max 



< a, 



where Or, 



a 



This results in a linear mass density 



-^^fil 



3a,; 



(3.22) 



for normal growth filaments.® The corresponding virial temperature is found in [60] to be 
Tvir ~ ^m-pG/ifii, where the proton mass nip = 1.1 x 10^^ K. Hence 



Gmp^mafi 



f-cq 



10 aj^_ 



a 
005 



0.3 



a,. 



K. 



(3.23) 
(3.24) 



Accelerated growth filaments will have a slightly smaller linear mass density. The difference will not be 
important except at very high redshifts when there are no star forming halos in the normal growth regime. 
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Thus star formation occurs only for loops formed after 

~ /x_8 V « J ^o-3>' yio^Kj 

Note that if T* = lO^K and the string tension obeys Gfi < 10~^, star formation can only be 
attributed to loops in the low-velocity tail of the loop production function, which is highly 

suppressed. 

If we restrict our attention to the ionization history of the universe, which is sensitive 
only to the total fraction of baryons in stars, we can neglect the subsequent dynamics of 
the filaments, which will fragment into beads which subsequently merge. This process will 
not significantly increase the virial temperature, since the filament is already in (2D) virial 
equilibrium, and the collapse into beads preserves the total energy. 

3.4 Longitudinal filament collapse (or lack thereof) 

Although we will find a linear instability of filaments toward collapse into bead-like halos, we 
can rule out the merging of the entire filament into a single large halo for all but the slowest 
loops. Here we consider the longitudinal collapse mode of the entire filament. 

To show such a filament will not collapse, we can disregard the finite size of the loop, 
since accounting for the finite size can only further delay collapse. Then within a Hubble 
time of teq the filament mass is roughly the mass of the loop, 

Mfii - afiti. (3.26) 



(3.25) 



The length of the filament at this time is Lfn ~ Veqteq, and it grows with the scale factor as 

Lm ^ Veqteqa, (3.27) 

where the loop velocity at teq is v^q ~ ctyai. The average mass per unit length is 

Mfii a/iOj 
/"fii ^ -J — ^ 

Let us first disregard the collapse along the filament axis. Then, at t > teq, both the 
mass of the filament and its length grow like (t/teq)^^^, so /Xfii remains constant. 

Now let us estimate the time scale of the longitudinal collapse. The characteristic 
longitudinal velocity that parts of the filament develop in a Hubble time t due to the filament's 
self-gravity is 



(3.28) 



V ~ (G/ifii/Lfii)i. (3.29) 

The collapse sets in when this becomes comparable to the Hubble velocity, vh ~ L/t. This 
gives the following estimate for the redshift of the collapse: 

(1 + Z) ~ aG^ia-^iteq/ti)^/\l + Zeq). (3.30) 

With Q: ~ 0.05, ay ~ 0.3, Gfi < 10"'^ and tj/tcq ^ 10~^, we get 1 + z < 1, i.e., no collapse. 
This means that the filaments from cosmic string loops remain elongated by a large factor. 



3.5 Fragmentation of filaments into beads 

As the filaments grow in mass and thickness from radial infall of dark matter (and baryons 
after recombination), the longitudinal expansion maintains their linear mass density at a 
constant (in time) value fJ-miu) = Peqxli^{y,a)/a. This overdense cylinder will be unstable to 
collapse into beads on scales longer than the Jeans length. The fastest growing instability 
was found in [62] to have a comoving wavelength Aj ~ 47rxvir, where the virial radius of 
a cylindrically symmetric self gravitating gas is Xvir = Xta/\/e ^ Xta, and so the comoving 
length of the typical bead is comparable to the transverse dimension of the cylinder, 

Ay ~ 4xta. (3.31) 

Longer wavelength instabilities represent merging of such beads, and occur on scales up to 
an order of magnitude longer. Since we characterize the beads by the size of the longest 
unstable mode, the length of the bead after merging is 



Lbead ~ WaAy ^ 207rtcqW « 4 x 10'^ y^a^a^fi^g Mpc. (3.32) 



The number of beads is then 



-^^bead 

In accelerated growth filaments, the bead mass is 



Z^eads — T"^- (3.33) 



M = ^ « 30a^/6/3-2/3^Meq ( , (3.34) 

i^cads \ Oiy J 

where Mgq = t^JG = 3.6 x lO^^Afo is (roughly) the mass contained in a horizon-sized sphere 

at teq* 

The mass of each halo from normal-growth filaments can likewise be estimated to be 
M=^^mMj^^"]'\ (3.35) 

By requiring T^i^. > T^,, we can set the lower bound on the relevant bead mass at time 
a{t) to be 

Mr - 50aiV6Meq./^ (^^Y' (3-36) 
\j nip \ ail3 ) 

for beads in accelerated growth filaments and 

/ rp \ 3/2 

^norm ^ I^M^^Ct'l'^ — (3.37) 



nip) 



for beads in normal growth filaments. 
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4 Halo spectrum 



We are interested in comparing models with cosmic strings of tension = 10~^/i-8 with 
the standard structure formation scenario. The halo spectrum of the standard scenario is 
well approximated by the Sheth-Tormen mass function [63]. The Acdm cosmology we use 
assumes the values of WMAP+BAO+Zfo [64], which for our purposes is just Zeq = 3232 and 
Mcq = 3.62 X IOI'^Mq. Because we are interested in early star formation, we will assume the 
growth function D{z) scales simply as D{z) tx (1 + 

We can calculate the spectrum of halos using the continuity equation. As we only 
consider radiation-era loops, there is no source term, since every loop already exists by 
ffl = Oeq = 1) and subsequently is associated with one filament, which in turn is associated 
with some number treads of halos. 

The spectrum of halos N{M, a) at any redshift z = ^"^^^ — 1 is determined entirely from 
the spectrum of loops n(a, m, v) at a = Oeq via 



N{a,M)= I l^hcads{a,m,Veq)n{aeq,m,Veq)-g^dVeq, 



where 



0.08v, 



7/6 fM^^^/^ /"G^V'^n-i/e 



2/3 



eq 



V P J 



< a < an 



< a, 



(4.1) 



(4.2) 



^eq 



a. 



2(2teq)3/2 (m + TGfiH 



cq_J 



x5/2 



and m(a; M) is given by combining (3.34) or (3.35) with the substitution 



a,; 



/ m 
\2afit, 



eq 



1/2 



(4.3) 



(4.4) 



The delta- function enables us to integrate (4.1), yielding V(,q — (m/2a/iteq)^/2. We can 
neglect the correction proportional to F, since the affected loops are too small to seed halos 
with Tvir > T*. 

The Jacobian {dm/dM) can be found from the loop mass 

11/3 



m{M) 



3.0 X 10 ^jjfa^ \a{Gp.)^a^ 



< a < On 



< a. 



(4.5) 



Note that although m{M) is continuous at the transition between the accelerated and normal 
growth regimes at o = amax, the derivative (dm/dM) is not, resulting in a discontinuity in 
the halo spectrum (4.1). 

We can write amin and amax from (3.19) without reference to by using (4.4), 

\ 3/2 



V2a( 

,2/5 



(4.6) 
(4.7) 
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We can also express amax in terms of the halo mass M by substituting the lower expression 
in (4.5) with a = amax into (4.6) and then solving for amax- This gives 



0.2a^ { (3M 



0.014 / M ^ ^''^ 



/i_8 V^G 



Gfi \a-Mcq 
The halo spectrum is then given by 

iVacc(a,M) a < ai 



(4.8) 



N{a,M) 



NnoT:m{a,M) a > Cmax, 



(4.9) 



where 



oa^BSV f Men 
2.5a5V f M 



a 



19/3 



iVnorm(a,M) 



Numerically, we find 



MiVacc(a,M)Mpcg = 4.9 x 10^*a^,/3(5P 



M7Vnorm(a,M)MpcE] = 2.8 x lO^^P— aV-g 



-5/3 



(4.10) 
(4.11) 

(4.12) 
(4.13) 



where the unit Mpco = Mpc/(1 + Zcq) is the comoving length which is a megaparsec today. 
We plot these distributions vs. Sheth-Tormen in figures 2 and 3 for three different values of 
Gfx. 

In the delta-function approximation for g, the largest contribution to the number density 
of halos is from the smallest loops (since their number density is larger). This is cut off when 
the resulting halos are too small to form stars. As we consider larger halos, we consider larger 
loops still in the normal growth regime. This is cut off when the loops are moving too fast 
to be in the normal growth regime. Even larger loops are responsible for accelerated growth 
halos. The small discontinuity in the graphs at the transition to this regime is caused by the 
discontinuity in the Jacobian factor {dm/ dM). The sharp decline of the halo distribution in 
the accelerated growth regime reflects the steep M~'^ dependence in (4.12). 

Prom the figures, there is a clear enhancement of early star formation provided G/i > 
lO^'' if T* = 10"^ K, and even for G^ > 10^'''^ if T* = 10^ K. This enhancement occurs for 
redshifts z > 10 — 20, and is limited to a rather narrow range of halo masses. Within the 
delta-function approximation for there is no early star formation for G/i < 10^^ ( io4*k ) ' 
although in fact, the low momentum tail of g will produce some early stars. 



5 Baryon collapse fraction 

The fraction of baryons which collapse into stars is proportional to the fraction of matter 
in star- forming halos, i.e., those of sufficient virial temperature for efficient atomic hydrogen 
cooling. Because the process of filament collapse does not significantly affect the virial tem- 
perature of the gas, we can neglect this process and simply count the fraction of matter in 
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5 10 15 4 6 8 10 12 14 



log,,, (M/Mg) log,o (M/Mg) 

Figure 2. Number and mass densities of lialos per logarithmic mass bin. Solid lines indicate halos 
with Tvir > lO'^K, and dotted lines extend this down to lO^K. The Sheth-Tormen mass function 
is used for comparison (thin lines). From right to left the redshifts are (1,6,10,20,50). We use the 
delta-function approximation for g with a = 0.05, a,j — 0.3, and dV — 7. 



filaments of sufficient mass to form stars, i.e., tliose of mass Mgi > Mgj. The collapse fraction 
is 



1 

Pm Jm, 
1 



HI 

oo 



(iMmMfiiiVm(Mfii), 
Mai(m)n(aeq, m, Veq)dmdveq, 



poo 
J nif 

SiraGSV y^ieqap^ 



2(2teq)3/2m5/2 

dm 



dm, 



m 



acc 



dm 

3/2 



5/2 



m 



dm 

5/2 



^acc ^ ^* 



(5.1) 
(5.2) 
(5.3) 

(5.4) 



where m=f, the smallest loop mass responsible for star forming filaments, is determined from 
Tvir > and (3.23) to be 



18MeqT2a2 



Gfiam"^ ' 

and the smallest loop mass undergoing accelerated growth is given by 

'2aa2(G;u)5^ 



(5.5) 



eq- 



(5.6) 
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5 10 15 4 6 8 10 12 14 

log,o (M/M^) log;,, (M/M^) 



Figure 3. Number density (left) and mass density (right) in halos with TVir > lO'^K (solid) and 
Tvir > lO^K (dotted) for strings of tension G^l = 10"'^ (top) and Gn = lO"^ '^ (bottom) . 



The filament mass Mai(m) from (3.20) then takes the form 

Mfii m = <^ - 5.7 

It will be convenient to define the time before which no star-forming filaments are in 
the normal-growth regime by m* = m^cc^ yielding 

_ MT^atP^ _ 1.4x IQS ^ T. ^ 

{G^ifmla^~ Mis VlO^Ky' • ^^'^^ 
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Then we can write 



Fool 



' ^TrGaSVy/ieqOfJ^ m-acc 

^ (5.9) 



It is here where we examine the consequences of using the delta-function approximation 
for the loop production function g. In particular, we have found that the large typical loop 
velocity does limit early star formation, and so we should verify that rare slow loops are not 
more important. Notice that the low-velocity tail of the distribution will rapidly diminish 
the value of Cacc oc a\. This means the normal-growth dominated (a > aacc) part of (5.9) is 
relevant, and so 

FeolOC<5P — . (5.10) 

Simulations reveal a loop production function g(a, Oy) consisting of several features, including 
a scaling peak and a broad plateau (see figures 5 &: 7 in [46]). This distribution becomes 
rather localized after multiplication by a and division by a^, since all features except the 
scaling peak become suppressed. The numerical results are well fit by a single localized peak 
for the distribution 

—g{a, ay) « 7^6{a - 0.05)<5(a^ - 0.3), (5.11) 

Oy 0.3 

and hence, for the purpose of calculating Fcoi, 

g{a, ay) ^ 7 6{a - 0.05)(5(a^ - 0.3). (5.12) 
For the potentially observable range of Gfj, and we can write 

'^aSV^iG,? ^ 6.5 X 10- (i°3) (5.13) 

2 ayl^ \ J 1 + z 



6 Ionization history 

The fraction of volume which has been re-ionized is related to the baryon collapse fraction 
Fcoi via [65, 66] 

where the number of ionizing photons per baryon in virialized halos A^ion depends on the 
mass distribution of stars. For the initial mass function of stars observed in nearby galaxies^, 
-^ion = 40(/^/0.1)(/esc/0.1), where is the efficiency of converting baryons in virialized halos 
into stars and /esc is the fraction of ionizing photons which escape the halo. Here 

C 

F(a',a) = -0.26— 
^ ' ^ 10 



^The value of A^ion increases by 1.5 orders of magnitude if the stars are all more massive than IOOMq [67]. 



where C = {n'^) / {nu)'^ is the clumpmess. We can evaluate Qhii(-z) in closed form using the 
incomplete Gamma function, although the result is rather long, so we will not present it 
explicitly here. The clumping factor C is expected to approach unity at high redshifts [68]; 
hence we set C = 1. We plot Fcoi{z) and Qmi{z) for various scenarios in figure 4 below 
usingS A^ion = 10 for = lO^K and iVion = 5 for = lO^K. 



^'coi(z) for T, = 10^ K F,M) for T, = lO'* K 




1 5 10 50 100 500 1000 \ 5 10 50 100 500 1000 

z z 



Figure 4. Fraction of matter in loop-seeded star-forming halos (top) and fraction of volume which 
has thereby been rc-ionizcd (bottom) as a function of redshift. We include a residual floor of fO^'' in 
Qhii- The Sheth-Tormen scenario is shown for comparison (thin black line). 

We can find the optical depth by integrating the reionization fraction out to redshift z 

via 

= nii{z')Quii{z')aT ^ z')H{z') ' ^^'^^ 

where n}i{z) is the physical number density of hydrogen, equal to 2.7 x 10~''cm^'^ today, 
and (Jt = 6.652 x 10~^^cm^ is the Thomson cross-section. The Hubble rate obeys H{z) = 
iJo\/f^A + f^m(l + zf, with JIa = 0.728 and = 1 - ^^A- 



These values of A'^ion are chosen for consistency with the WMAP optical depth; sec figure 7 below. 
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The visibility function V{z) is defined through the optical depth r via, 

Viz) = e--^^^H{z)^. (6.4) 
We plot t{z) and V{z) in figures 5 & 6. 



t(z) for = 10^ K t(z) for T, = 10'' K 




Figure 5. Optical depth from Shcth-Tormen plus cosmic strings for two possible values for T*. 



V{z) Mpc for = 10^ K V{z) Mpc for = lO"* K 




Figure 6. Visibility function from Sheth-Tormen plus cosmic strings for two possible values for T*. 

Using the WMAP+BAO+ffo [64] measured value of the reionization optical depth 
r = 0.087 ± 0.014, we can find the excluded region^ of Nion-Gfj, parameter space, shown in 
figure 7 below. 

Using the functional form of Qmiiz), we examine the effects of early reionization due to 
strings on the CMB temperature and ii^-mode polarization power spectra in figures 8 & 9, as 
well as on the TE cross-correlation in figure 10. In these figures we plot both the pure Sheth- 
Tormen case as well as Sheth-Tormen plus cosmic strings for the specified values of Gfi and T*. 
(Note that here we have accounted only for the effect of strings on reionization, disregarding 

®The parameters indicated in figure 7 are excluded only within the framework of our simplified model, 
which assumes that A^ion, T,, and C are constants independent of the redshift. In particular, larger values of 
A'^ion niay be allowed at smaller redshifts, where the dumpiness C is likely to be > 1. 
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their direct gravitational influence on CMB spectra, e.g., through the Gott-Kaiser-Stebbins 
effect.) We additionally plot the fractional difference in power spectra induced by cosmic 
strings. These figures were obtained using a modification of the CAMB^*^ software package 
[69]. Because Planck is expected to be near cosmic variance limited [70], we have shaded 
the regions of these figures below the cosmic variance detection threshold for a single Ci. 
(Features consisting of several consecutive Qs can be estimated more accurately.) 



Contour for WMAP+BAO+Z/p: r = 0.087 + 0.014 



Contour for WMAP+BAO+i/n: r = 0.087 ± 0.014 





5 10 15 20 5 10 15 20 25 30 

10** Gyu 10** Gn 

Figure 7. Predicted optical depth vs. A^ion and G/i. Parameter space shaded with gray is compatible 
with the 68% CL WMAP+BAO+i?o [64] reionization optical depth. Note that our value for optical 
depth is defined by integrating from z = to ^decoupling- The pure Shcth-Tormen scenario is equivalent 
to Gn = 0. 



7 Discussion 

Loops of cosmic string can seed relatively massive halos at very early times. These halos can 
become sites of early star formation and can influence the ionization history of the universe. 
Figures 8-10 illustrate that for Gfi > 10^^ the resulting effect on the CMB temperature and 
polarization power spectra can be significant and is likely to be detectable by Planck. Given 
the present uncertainties in the physics of reionization, we did not attempt to translate our 
analysis into a firm bound on G/i. 

As we mentioned in the introduction, values of G/i > 10"'' are in conflict with the bound 
G/i ^ 4 X 10~^, obtained from millisecond pulsar measurements in Refs. [25, 27]. This bound, 
however, does not account for the fact that a significant fraction of the string network energy 
goes into the kinetic energy of loops, which is then redshifted. It also assumes that nearly all 
loops develop cusps, while examination of loops produced in the recent simulations suggests 
that cusps are rather rare. Accounting for these effects is likely to relax the pulsar bound by 
at least an order of magnitude. 

1° http://camb.info/ 
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Gn^O and Gft = 10"' for T, = IO'K G// = and Gf^ = lO"**"' for T, = lO'^K 




Figure 8. CMB temperature angular power spectra with and without cosmic strings (top) and their 
fractional difference (bottom). We assume T* — lO^K (left) and T* = lO^K (right). The shaded 
region indicates the single C/ cosmic variance for the G/z — case. 



We note that the current census of observed star formation at redshifts 6-10 falls short 
of accounting for the minimum production rate of UV photons that is required for ionizing 
the universe [71, 72]. A string-seeded scenario might be called for if future observations would 
reveal an abundant population of galaxies at redshifts z > 10. In the standard cosmological 
model for structure formation, the comoving star formation rate is expected to decline sharply 
at increasing redshifts beyond z ~ 10, since the massive halos in which gas can cool and 
fragment into stars become progressively rarer on the exponential tail of the Gaussian density 
field. The cosmic string scenario provides a natural way to account for additional star 
formation at high redshifts, z > 10. 

Early star formation from cosmic strings may help explain the observed metallicity of 
the intergalactic medium (IGM) at redshifts z > 2. In particular, the observed metallicity 
floor of C IV at redshifts as high as z ~ 6 [73] cannot be explained with the standard star 
formation history in which the filling fraction of enriched IGM is very small at these redshifts 
[74]. 

Additional implications of early star formation in the cosmic string scenario involve 
the existence of an abundant population of gamma-ray bursts at unusually high redshifts, 

^^Currently, the earliest observed GRB is at z ~ 9.4 [75], which is within the range expected from the 
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G// = and Gfi = 10"^ for T, = lO^K 



Gyt/ = and Gfi = 10"''^ for T, = IO^'K 




Figure 9. CMB i?-mode polarization angular power spectra with and without cosmic strings (top) 
and their fractional difference (bottom). We assume T^. = lO^K (left) and T^, = 10*K (right). The 
shaded region indicates the single C; cosmic variance for the G/x = case. 



as well as early formation of massive galaxies, which could be detected by future infrared 
telescopes such as JWST^^, and early supermassive black-holes^'^, which could be detected 
by future X-ray telescopes such as 1X0/ Athena^'* or gravitational wave observatories such as 

Early structure formation around cosmic string loops could also be tested by future 
21-cm observations [76]. For example, "global" 21-cm measurements, which give the sky- 
averaged spectrum of the 21-cm signal, are sensitive to early production of Lya photons 
(which couple the 21-cm excitation temperature of hydrogen to its kinetic temperature) at 
z < 50 and to the subsequent production of X-rays which heat the intergalactic hydrogen. 
These ingredients combine to make a distinctive absorption trough in the global 21-cm signal, 
which could potentially be measured by future experiments [77]. 

In this paper we focused on "ordinary" cosmic strings, characterized by Nambu-Goto 

standard cosmological model, 
^^http:/ /www .jwst.nasa.gov/ 

^^The possibility of early super-Eddington growth of black holes in high redshift galaxies is reviewed in 
Ref. [78]. 

^"^http:/ /sci.esa.int/ixo 
^^http:/ /sci.esa.int/lisa 
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Gyu = and 10"^ for = lO^K G// = and lO-*"^ for T, = IO^'K 
2.0, . . . 1 2.0| . . . 




Figure 10. CMB TE cross-correlation with and without cosmic strings (top) and their relative 
difference (bottom). We assume = fO^K (left) and T* = lO'^K (right) for the indicated tensions 
vs. pure Sheth-Tormen. The shaded region between ±^2/{2l + 1) represents the Gfi = cosmic 
variance limit for a single C;. 

equations of motion and reconnection probability p^ec = 1 • The value of prec could be much 
smaller for cosmic superstrings, resulting in a denser string network, a higher density of 
loops, and a larger effect on structure formation (for the same value of G/x). Superstrings 
can also form Y-junctions, with three strings meeting at a vertex. The evolution of networks 
with junctions and/or with prec < 1 is quantitatively different from that of ordinary strings, 
and reliable observational predictions would require large-scale numerical simulations of such 
networks. 
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A The rocket effect 



The gravitational radiation from a cosmic string loop is anisotropic in general. This affects 
the momentum of the loop. Let the radiation be beamed in the minus y-direction. The 
physical y-momentum p = mv / \/\ — of the loop then obeys 

| = -Fp + rpG/i2, (A.i) 

whereas Tp ~ 10. 

If we assume the non-relativistic loop has velocity Ueq > in the +2/-direction at teq, its 
subsequent velocity is given by 

_ ^ ^ 3rp(Grt^M,,(„V^-a-')_ 

a 5m 

where a = (f/teq)^^^- The trajectory in comoving coordinates is then 

9rp(GA^)^Meqteq ^ , 2 ^ o , ^ M f ^ 

y{a) = — ^-^ + a - 5j + 3.eqieq " ^ J ' (A.3) 

and so the matter-era velocity of the loop formed at in the radiation era is 

v{a, a,) = Q.- + ^ 2 ~ a.- + ^ (^-4) 

From this, we can calculate the turnaround surfaces, and find the average filament 
thickness by 

which determines the average bead mass found in section 3.5. For filaments from loops 
produced after Oj ~ 0.02^_8, this leads to a correction factor of order 0.8 to the average 
filament radius. The only loops whose filaments receive significant corrections from the rocket 
effect are those whose number-densities are highly suppressed due to loop evaporation. Loop 
evaporation becomes significant only for loops of mass below 

mp = r(G^)2Meq. (A.6) 

Since no star formation results from loops of mass less than given by (5.5), we can neglect 
loop evaporation for 

* < 1 (AT) 

25{G^ifamj- ' ^ 

i.e., 

/ rp \ 2/3 

For this reason, we can neglect loop evaporation and the rocket effect entirely. 



^®This value was estimated for simple Kibble- Turok [79] loop solutions involving a few lowest harmonics. 
Recent (non-scaling) simulations have shown that loops chopped off the network tend to be rectangular [80], 
which can suppress Fp. 
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